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Abstract 

A search for Higgs boson decays to invisible particles is performed using 20.3 fb”^ of pp 
collision data at a centre-of-mass energy of 8 TeV recorded by the ATLAS detector at the 
Large Hadron Collider. The process considered is Higgs boson production in association 
with a vector boson (V - W or Z) that decays hadronically, resulting in events with two or 
more jets and large missing transverse momentum. No excess of candidates is observed in 
the data over the background expectation. The results are used to constrain VH production 
followed by H decaying to invisible particles for the Higgs boson mass range 115 < mn < 
300 GeV. The 95% confidence-level observed upper limit on (Tvh x BR{H —> inv.) varies 
from 1.6 pb at 115 GeV to 0.13 pb at 300 GeV. Assuming Standard Model production and 
including the gg H contribution as signal, the results also lead to an observed upper limit 
of 78% at 95% confidence level on fhe branching rafio of Higgs bosons decays fo invisible 
particles at a mass of 125 GeV. 
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1 Introduction 


Since the discovery of a Higgs boson with a mass of approximately 125 GeV [1, 2] at the LHC in 2012, 
the properties of this new particle have been studied extensively. All results obtained so far [3-9] are 
consistent with the expectations of the long-sought Standard Model (SM) Higgs boson [10-13]. However, 
sizeable deviations from the SM expectation cannot be yet excluded; the total branching ratio of beyond- 
the-SM decays of the Higgs boson is only weakly constrained, and its value could be as high as ~ 40% [8, 
14]. One possible decay is to weakly interacting particles, as predicted by many extensions of the SM, e.g. 
Higgs boson portal models [15-18]. In these models, the Higgs boson can decay to a pair of dark-matter 
particles if kinematically allowed. These decays are generally “invisible” to detectors, resulting in events 
with large missing transverse momentum (L"™®®). 

Searches for Higgs boson decays to invisible particles {H inv.) have been performed by both the 
ATLAS and CMS collaborations [14, 19]. For example, the ATLAS Collaboration has placed an upper 
limit of 75% [19] on the branching ratio of // —> inv. from Higgs boson production in association with 
a Z boson identified from its leptonic decays (Z —> The present paper describes an independent 

search for the H inv. decay in final sfafes wifh fwo or more jefs and large F™®®, mofivafed by Higgs 
boson producfion in association wifh a vector boson V {V = IT or Z): qq' —> VH. The vector boson 
is idenfified fhrough ifs decay fo a pair of quarks, reconsfrucfed as hadronic jefs in fhe ATLAS defecfor, 
V jj. Gluon fusion production gg ^ H followed by // —> inv. can also lead fo evenfs wifh fwo or 
more jefs and large F™*®, and Iherefore confribufes to fhe signal of fhe search. Negligible confribufions 
of approximafely 1% and 0.2% to fhe sensifivify come from qq' —> qq'H producfion via vecfor-boson 
fusion (VBF) and from qqigg —> ttH (ttH) producfion, respectively. The VBF confribufion is sfrongly 
suppressed by fhe mjj (dijef invarianf mass) window cufs and by fhe forward-jel veto used fo reduce fhe lop 
quark-anliquark background (ti), as described in Seel. 4. In a previous ATLAS dark-malfer search, limils 
on Higgs boson decays fo invisible particles in VH production were sel using evenfs wifh a hadronically 
decaying veefor boson and F™®® as well [20]. However, fhe presenl analysis achieves better sensifivify by 
using differenl fechniques and performing dedicaled oplimizalions. 


2 Experimental setup 


This search is based on prolon-prolon collision dala al a cenlre-of-mass energy of 8 TeV recorded wifh 
fhe ATLAS defector [21] in 2012, corresponding to an integrated luminosily of 20.3 fb“^ The AT¬ 
LAS defecfor is a general-purpose defector wifh an inner fracking syslem, eleclromagnefic and hadronic 
calorimefers, and a muon speclromeler surrounding fhe inleraclion poinf. * The inner fracking system is 
immersed in a 2 T axial magnefic field, and fhe muon speclromeler employs a toroidal magnelic field. 
Only dafa recorded when all subdelecfor sysfems were funcfional are used in Ihis analysis. 

The Irigger syslem is organised in Ihree levels. The firsl level is based on custom-made hardware and 
uses coarse-granularify calorimeter and muon informalion. The second and Ihird levels are implemented 
as soflware algorilhms and use fhe full defector granularify. Al fhe second level, only regions deemed 


* The ATLAS experiment uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the 
centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the 
i/-axis points upward. Cylindrical coordinates (r, (p) are used in the transverse plane, 0 being the azimuthal angle around the 
z-axis. The pseudorapidity is defined in terms of the polar angle 0 as 77 = - ln[tan(6/2)]. Transverse momenta are computed 
from the three-momenta, p, as pj = |/;| sin S. 
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interesting at the first level are analysed, while the third level, called the event filter, makes use of the full 
detector read-out to reconstruct and select events, which are then logged for offline analysis at a rate of 
up to 400 Hz averaged over an accelerator fill. 


3 Object reconstruction and simulated samples 

Jets are reconstructed using the anti-kf algorithm [22] with a radius parameter of = 0.4. Jet energies are 
corrected for the average contributions from minimum-bias interactions within the same bunch crossing 
as the hard-scattering process and within neighbouring bunch crossings (pile-up). Furthermore, for jets 
with pt < 50 GeV and \r]\ < 2.4, the scalar sum of the pj of tracks matched to the jet and originating 
from the primary vertex^ must be at least 50% of the scalar sum of the pj of all tracks matched to the 
jet, to suppress jets from pile-up interactions. Jets must have pj > 20 GeV {pj > 30 GeV) for I 77 I < 2.5 
(2.5 < |77| < 4.5). 

Jets containing fj-hadrons (ft-jets) are identified (Z^-tagged) using the MVlc algorithm, which is an im¬ 
proved version of the MVl algorithm [23] with higher rejection of jets containing c-hadrons (c-jets). It 
combines in a neural network the information from various algorithms based on track impact-parameter 
significance or explicit reconstruction of secondary decay vertices. The operating point of this algorithm 
chosen for this analysis has an efficiency of about 70% for f^-jets in tt events and a c-jet (light-jet) mis-tag 
rate less than 20 % ( 1 %). 

Lepton (electron or muon) candidates are identified in two categories: loose and tight, in order of increas¬ 
ing purity. Electron candidates are reconstructed from energy clusters in the electromagnetic calorimeter 
matched to reconstructed tracks in the inner tracking system. They are identified using likelihood-based 
methods [24, 25] . Loose electrons must satisfy “very loose likelihood” identification criteria and are re¬ 
quired to have pi > 1 GeV and \r]\ < 2.47. Tight electrons are selected from the loose electrons and 
must also satisfy the “very tight likelihood” identification criteria. Muon candidates are reconstructed 
using information from the inner tracker and the muon spectrometer [26]. Loose muons are required to 
have Pj > 7 GeV and \t]\ < 2.7. Tight muons are then selected from the loose muons, by requiring 
Pj > 25 GeV and \r]\ < 2.5. They must be reconstructed in both the muon spectrometer and the inner 
tracker. Lor the loose leptons, the scalar sum of the transverse momenta of tracks within a cone of size 
AR = + {AtjY = 0.2 around the lepton candidate, excluding its own track, is required to be less 

than 10% of the transverse momentum of the lepton. Lor the tight leptons, there are more stringent isol¬ 
ation requirements: the sum of the calorimeter energy deposits in a cone of size A/? = 0.3 around the 
lepton candidate, excluding the energy associated with it, must be less than 4% of the lepton candidate 
energy, and the track-based isolation requirement is tightened from 10% to 4%. 

The missing transverse momentum vector, is computed using fully calibrated and reconstructed 

physics objects, as well as clusters of calorimeter-cell energy deposits that are not associated with any 
object [27]. Only calibrated jets with pj greater than 20 GeV are used in the computation. The jet energy 
is also corrected for pile-up effects [28]. A track-based missing transverse momentum vector, is 

calculated as the negative vector sum of transverse momenta of reconstructed tracks associated with the 
primary vertex and within I 77 I < 2.5. 

Monte Carlo (MC) simulated samples are produced for both the signal and background processes. Unless 
otherwise stated, the simulation [29] is performed using the ATLLAST-II package [30], which combines 

^ The primary vertex is taken to be the reconstructed vertex with the highest of the associated tracks. 
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Table 1: List of MC generators, parton distribution functions (PDFs) and cross sections used for the signal and 
background processes. The H —» inv. signal cross sections are given for nin - 125 GeV and assume SM production 
and BR(// —> inv.) = 100%. Details are given in the text. 


Process 

Generator 

PDEs 

Cross section [pb] 

tt 

Powheg -i- Pythia 

CTIO [42] 

normalized to data 

F-bjets 

Sherpa 

CTIO 

normalized to data 

Single top 

?-channel 

AcerMC 

CTEQ6E1 [43] 

88 

5-channel 

Powheg -i- Pythia 

CTIO 

5.6 

Wt 

Powheg -i- Pythia 

CTIO 

22 

Diboson 

ITIT 

Powheg -i- Pythia 

CTIO 

52 

WZ 

Powheg -i- Pythia 

CTIO 

9.2 

ZZ 

Powheg -i- Pythia 

CTIO 

3.3 

SM VH 

qq' VH{-^ bb) 

Pythia 

CTEQ6E1 

0.18 

gg ZH{-^ bb) 

Powheg -i- Pythia 

CTIO 

0.0038 

Signals 

qq Z{-^ jj)H{-^ inv.) 

Herwig-p-i- 

CTIO 

0.29 

qq' ^ jj)H{^ inv.) 

Herwig-p-i- 

CTIO 

0.48 

gg H{-^ inv.) 

Powheg -p Pythia 

CTIO 

19 


a parameterized simulation of the ATLAS calorimeter with the GEANT4-based [31] full simulation for the 
rest of the subdetector systems. 

Signal events from qq' —> VH with H inv. are produced using the NLO Powheg method as implemen¬ 
ted in the Herwig-i-i- generator [32]. The gg ZH production process contributes approximately 5% 
to the total ZH cross section. Events from the gg ZH production process are not simulated, but are 
taken into account by increasing the qq ZH cross section as a function of the Higgs boson pj by the 
appropriate amount. The gluon-fusion signal events are produced using the Powheg generator interfaced 
to PythiaS for parton showering and hadronization. The production of qq' —> VH followed by the SM 
H ^ bb decay is considered as a background for the search. The PythiaS generator is used to produce 
these events. The cross sections of all Higgs production processes are taken from Ref. [33]. 

A significant source of background is the production of P-i-jets and of tt events. A sample of P-i-jets events 
is generated using the Sherpa generator [34] with massive b- and c-quarks. Events from the tt process 
are generated using the Powheg generator interfaced with Pythia6 [35]. Other background contributions 
include diboson (WW, WZ and ZZ) and single top-quark production. The Powheg generator interfaced to 
PythiaS is used to produce diboson events. The diboson cross sections are calculated at NEO in QCD 
using the MCEM program [36] with the MSTW2008NEO parton distribution functions (PDEs) [37]. The 
5 -channel and Wt single top-quark events are produced using the Powheg generator, as for tt production. 
The remaining t-channel process is simulated with the AcerMC generator [38] interfaced to Pythia6. 
Cross sections of the three single top-quark processes are taken from Refs. [39-41]. Table 1 summarizes 
the MC generators, PDEs and normalization cross sections used in this analysis. 
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4 Event selection 


Events are required to pass an trigger with a threshold of 80 GeV, which is a cut applied at the third 
level. The E™*® trigger is fully efficient for E™*® >160 GeV and 97% efficient for E™®® = 120 GeV. An 
efficiency correction is derived from W //v+jets and Z —> //"'■//“+jets events. This correction is below 
1% for 120 GeV< E™®® <160 GeV. Events are also required to have >120 GeV, p™*® > 30 GeV, 
no loose leptons and two or three “signal jets” (satisfying \r]\ < 2.5, px > 20 GeV and leading jet 
Pt > 45 GeV). The inclusion of 3-jet events improves the signal efficiency. A requirement is made on 
Hj, defined as the scalar sum of the px of all jets: Hj > 120 (150) GeV for events with two (three) 
jets. This cut is employed to avoid a trigger bias introduced by the dependence of the trigger efficiency 
on the jet activity, as also discussed in Ref. [44]. Events are discarded if they have additional jets with 
Px > 20 (30) GeV and \r]\ < 2.5 (2.5 < \t]\ < 4.5) to reduce the contribution from the tt background 
process. 

Eor VH signal events, E™*® resulting from the H —> inv. decay is expected to be strongly correlated with 
the transverse momentum of the vector boson V (pj). Since the E™** distribution of the signal is harder 
than that of the background, additional sensitivity in the analysis is gained by optimizing the selection 
cuts separately for four E™®* ranges. Here and in the following, the dijet refers to the two leading jets 
in events with three jets. The dijet invariant mass, mjj, is required to be consistent with that of the W/Z 
boson. In addition a requirement on the radial separation between the two jets, AEp, is made as the jets 
are expected to be close in for highly boosted V-bosons. Both the mjj and the AE// cuts reduce the V-i-jets 
and the tt backgrounds, and depend on E™'®*. The cut values are given in Table 2. 

Multijet events are copiously produced in hadron collisions. Eluctuations in jet energy measurements 
in the calorimeters can create E™'®* in these events and therefore mimic the signal. To suppress their 
contribution, additional selection criteria are applied to the azimuthal angles between p™*® and 

jets: A0(E™^^p™*^) < n/l, min[A0(E™^^jet)] > 1.5 and A0(E™^^ dijet) > 2.8. Here A0(E™^^p™^^) is 
the azimuthal angle between E™'®* and p™®®, min[A0(E™*®, jet)] the angle between E™'^® and its nearest 
jet, and A0(£™®®, dijet) is the angle between £™®® and the momentum vector of the dijet system. These 
requirements are based on characteristics of events with mismeasured E™'^® in the multijet background, 
while taking advantage of the expected topologies of signal events. 

Einally, the selected events are further categorized according to b-tag multiplicity (zero, one and two 
(7-tagged jets) to improve the sensitivity. Combined with the two categories in jet multiplicity (two and 
three jets), there are in total six categories in the signal region. 


Table 2: The E™***-dependent event selections of the signal region for the four E™'*’* ranges. 


E“‘®® range [GeV] 

120-160 

160-200 

200-300 

>300 

Variable 


Selection 


AEp, 2- and 3-jet events 

0.7-2.0 

0.7-1.5 

< 1.0 

<0.9 

Mjj, 2-jet events [GeV] 

70-100 

70-100 

70-100 

75-100 

ntjj, 3-jet events [GeV] 

50-100 

55-100 

60-100 

70-100 
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5 Background estimation 


In addition to the signal region, a number of control regions, designed to estimate various background 
contributions, are defined. They include the signal sideband (events not passing the mjj requirement), and 
the regions dominated by T+jets and tt events as discussed below. The multijet background is estimated 
from the data. The distributions of the T+jets and ft backgrounds are taken from MC simulation while 
their normalizations are estimated from the data. The remaining diboson, single-top and SM VH(bb) 
backgrounds are obtained from MC simulation. 

The multijet background is estimated using four regions defined by requirements on and 

min[A0(£'™*®, jet)], as listed in Table 3. The shapes of the rujj and distributions in the signal region 
A are taken from region C and the normalizations are determined by the ratio of the numbers of events in 
regions B and D. 


Table 3; Definition of the signal region, A, and the three regions B, C and D used to estimate the multijet background 
in the signal region. 


Region 

A 

B 

C 

D 


< n/2 

< n/2 

> n/2 

> n/2 

min[A0(£™*^ jet)] 

> 1.5 

< 0.4 

> 1.5 

<0.4 


The normalizations of the T-i-jets backgrounds are estimated using control regions enhanced in IT-i-jets 
and Z-rjets events. In all cases at least one lepton is required to have pj >25 GeV. The IT-i-jets events 
are selected by requiring exactly one tight lepton, > 20 GeV (S™®® > 50 GeV if Pj > 200 GeV), 
exactly two signal jets and < 120 GeV.^ Moreover, pj > 100 GeV is required in order to approxim¬ 
ately match the phase space of the signal region. The Z-i-jets events are selected by requiring two loose 
leptons of the same flavour with opposite charges with invariant mass 83 < nia < 99 GeV, at least two 
signal jets and a dilepton transverse momentum greater than 100 GeV. The kinematic distributions of the 
V-i-jets backgrounds are obtained from simulation that takes into account the different flavour composition 
of the jets. The simulated events are reweighted depending on the A^(jetj, jet 2 ) and to better match 
the data distributions [44]. The Z-rjets control region has a small contribution from ti (1.3%), which is 
estimated using a tt control region. This region is selected by requiring events to have two oppositely 
charged leptons of different flavour (one of which has pj >25 GeV) and passing the loose selection re¬ 
quirements, and at least two signal jets which are Z^-tagged. The signal sideband and the V-i-jets control 
regions are divided to match the categorization of the signal region while the tt control region remains as 
one category as described above. For the V-i-jets and tt control regions, the distributions of the multijet 
background are obtained from control regions defined by inverting the lepton isolation requirement and 
the normalizations are determined by template fits [44] . 


^ The transverse mass, m^, is calculated from the transverse momentum and the azimuthal angle of the charged 
lepton, Pj and and from the missing transverse momentum’s magnitude, £'™“, and azimuthal angle, = 

— cos((f>^ — 0™*®)). The transverse momentum of the W boson, Pj , is reconstructed as the magnitude of the 
vector sum of the lepton transverse momentum and the 
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6 Systematic uncertainties 


The experimental systematic uncertainties considered include the trigger efficiency, object reconstruction 
and identification efficiency, and object energy and momentum scales as well as resolutions. Among 
these, the jet energy scale (JES) and resolution (JER) uncertainties have the largest impact on the result. 
The JES uncertainties are ±3% and ±1% for central jets with a pj of 20 GeV and 1 TeV, respectively. The 
JER uncertainty varies from between ±10% and ±20%, depending on the pseudorapidities of the jets, for 
jets with pt = 20 GeV to less than ±5% for jets with px > 200 GeV. The JER and JES uncertainties are 
also propagated to the E™®* uncertainty. The f?-tagging uncertainty depends on jet px and comes mainly 
from the uncertainty on the measurement of the efficiency in tt events [23]. The dominant contribution 
arises from jets matched to fj-hadrons in the MC record of the particles’ true identities. Their efficiency 
uncertainties are at the level of ±2-3% over most of the jet px range, but reach ±5% for px = 20 GeV 
and ±8% above px = 200 GeV [45]. The uncertainty on the integrated luminosity is ±2.8%. It is derived 
following the same methodology as that detailed in Ref. [46]. 

Eor the backgrounds, a large number of modelling systematic uncertainties are considered, which account 
for possible differences befween fhe dafa and fhe MC models. These uncerfainfies are estimated following 
fhe sfudies of Ref. [44] and are briefly summarized here. The uncerfainfies on fhe V±jefs backgrounds 
come mainly from fhe knowledge of jef flavour composifion and fhe p^, and nijj disfribufions. Eor tt 
producfion, uncerfainfies on fhe fop quark fransverse momenfum and fhe nijj, E™'*® and p^ disfribufions 
are considered. The diboson background uncerfainfies are dominated by the theoretical uncertainties 
of the cross-section predictions, which include contributions from the renormalization and factorization 
scales and the choice of PDEs. The robustness of the multijet background estimation is assessed by 
varying the definition of the control regions B and D and an uncertainty of ±100% is assigned for this 
small background (< 1% in the signal regions). 

The uncertainty on the signal acceptance is evaluated by changing the factorization and renormalization 
scale parameters, parton distribution function choices and the parton shower choices. Eor the VH signal, 
the dominant uncertainty is from parton shower modelling, which can be as large as ±8%. Eor the gg ^ H 
signal, the dominant uncertainty originates from the renormalization and factorization scales and can be 
as large as ±15% in the high E™®® regions. Additional corrections to the Higgs boson px distribution of 
the gg ^ H signal are applied to match the distribution from a calculation at NNLO±NNEE provided by 
HRes2.1 [47, 48]. The detailed precedures are following the ones used in the H ^ yy and H —> WW* 
analyses as described in Ref. [49, 50]. The related uncertainties are also taken into account. 


7 Results 

The potential H inv. signal is extracted through a combined likelihood fit to the observed E™®® distri¬ 
butions of the signal region and its sideband and the p^ distributions of the control regions (p^ is defined 
as p^, p| and p^^^ for fhe lT±jefs, Z±jefs and tt confrol regions, respectively). The normalizafions of 
fhe V±jels and tt backgrounds are free parameters in this fit. The E™®* distributions are binned in such a 
way that each bin yields approximately the same amount of expected signal. The 2-jet categories of the 
signal region are split into ten bins, while fewer bins are used in the 3-jet categories and the sideband. 
Most V±jets control regions are split into five p^ bins, each yielding approximately the same amount of 
expected background. The 0-tag category of the V±jets control regions and the tt control region are used 
inclusively in the fit. The signal strength p, defined as the ratio of the signal yield (cry// x BR(// ^ inv.)) 
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Table 4: Predicted and observed numbers of events for the six categories in the signal region. The yields and 
uncertainties of the backgrounds are shown after the profile likelihood fit to the data. In this fit all categories share 
the same signal-strength parameter. The quoted uncertainties combine the statistical and systematic contributions. 
These can be smaller for the total background than for individual components due to anti-correlations. The yields 
and uncertainties of the signals are shown as expected before the fit for niff = 125 GeV and BR(H —> inv.) = 100%. 


Signal contributions from VBF and ttH production are estimated to be negligible. 


b-tug category 

0-tag 


1-tag 


2-tag 



Process 




2-jet events 




Background 

Z-i-jets 

24400 

+ 

1100 

1960 

+ 

200 

164 

+ 

13 


VT-i-jets 

20900 

+ 

770 

1160 

+ 

130 

47 

+ 

7 


tt 

403 

+ 

74 

343 

+ 

65 

57 

+ 

10 


Single top 

149 

+ 

16 

107 

+ 

14 

11 

+ 

2 


Diboson 

1670 

+ 

180 

227 

+ 

25 

64 

+ 

7 


SM NY^ibb) 

1.5 

+ 

0.5 

6 

+ 

2 

3 

+ 

1 


Multijet 

26 

+ 

43 

8 

+ 

7 

0.7 

+ 

0.9 

Total 


47560 

+ 

490 

3804 

+ 

64 

347 

+ 

15 

Signal 

gg ^ H 

403 

+ 

95 

25 

+ 

6 

2.1 

+ 

0.5 


iiW 

425 

+ 

45 

44 

+ 

6 

0.6 

+ 

0.1 


ji)H 

217 

+ 

19 

42 

+ 

4 

26 

+ 

2 

Data 


47404 


3831 



344 


3-jet events 

Background 

Z-i-jets 

9610 

+ 

580 

795 

+ 

93 

53 

+ 

7 


VT-i-jets 

7940 

+ 

510 

479 

+ 

70 

21 

+ 

4 


tt 

443 

+ 

53 

437 

+ 

53 

63 

+ 

7 


Single top 

97 

+ 

14 

66 

+ 

9 

6.4 

+ 

0.9 


Diboson 

473 

+ 

54 

55 

+ 

6 

13 

+ 

2 


SM NYi{bb) 

0.8 

+ 

0.3 

2.6 

+ 

0.9 

1.4 

+ 

0.5 


Multijet 

22 

+ 

29 

4 

+ 

4 

0.6 

+ 

0.6 

Total 


18580 

+ 

200 

1840 

+ 

40 

158 

+ 

7 

Signal 

gg ^ H 

224 

+ 

55 

15 

+ 

4 

1.2 

+ 

0.5 


iiW 

110 

+ 

16 

11 

+ 

1 

0.14 

±0.03 


ii)H 

65 

+ 

7 

12 

+ 

1 

6.1 

+ 

0.7 

Data 


18442 


1842 



159 




relative to the SM production cross section and assuming BR{H inv.) = 100%, is used to parameterize 
the signal in the data. A binned likelihood function is constructed as the product of Poisson probability 
terms comparing the numbers of events observed in the data to those expected from the assumed signals 
and estimated background contributions for all categories of the signal and control regions. The likeli¬ 
hood takes into account the background normalization and the systematic uncertainties. It is maximized 
to extract the most probable signal-strength value, ju. 

Table 4 shows the numbers of observed events in the data compared to the numbers of estimated back¬ 
ground events from the likelihood fit for each signal category. In all categories the data agrees with the 
background estimation. The backgrounds are dominated by Z-i-jets and IT-f-jets events. Subleading back¬ 
grounds come from top and diboson production. The SM VH and multijet background contributions are 
very small with the final even! selecfion. 

The fif reveals no significanf excess of evenfs over fhe background expecfafions and yields a besl-fil signal- 
slrengfh value of fi = which is consisfenf wifh zero. The confribufions from fhe individual 

sysfemafic uncerfainfies are summarized in Table 5. The sysfemafic uncertainly sources which have 
fhe largesl impacls are fhe energy scale of fhe jels and of along wifh fhe modelling (shape and 

normalization) of fhe diboson and P-i-jels backgrounds. The disfribulions of fhe evenfs passing fhe 

signal region selecfion are shown in Figs. 1 and 2 afler fhe profile likelihood fil fo fhe dafa. The fil resulfs 
are also propagafed fo fhe mjj disfribulions of fhe evenfs passing fhe signal region selection (wilhoul fhe 
m^y-window culs). The corresponding plols are shown in Figs. 3, 4 and 5 for fhe fhree b-tag cafegories 
separalely. 


Table 5: Impacts of sources of systematic uncertainty on the uncertainty of the fitted signal strength, A/j., in the 
data. Only sources with contributions larger than +0.03 are listed. 


Source 

Impacl 

on Aji 

Objecl sysfemafic uncerfainfies 

Jefs & 

+0.22 

-0.22 

Luminosily 

+0.04 

-0.03 

(^-lagging 

+0.05 

-0.04 

Background sysfemafic uncertainties 

Diboson 

+0.26 

-0.29 

Z+jels 

+0.21 

-0.22 

VF+jels 

+0.15 

-0.16 

tt 

+0.06 

-0.05 

Mulfijel 

+0.07 

-0.07 

Tolal 

Tolal sysfemafic uncertainly 

+0.41 

-0.43 

Dafa slalislical uncertainly 

+0.12 

-0.12 

Tolal uncertainly 

+0.43 

-0.44 


The null results are used to set 95% confidence level (CL) upper limils on fhe producl of fhe VH cross 
seclions and fhe V jj and H inv. decay branching rafio, (Tvh x BR(// —> inv.), as a function of fhe 
Higgs boson mass in fhe range 115 < mn < 300 GeV as shown in Fig. 6. The limits are computed with 
a modified frequenfisl melhod, also known as CL^ [51], and a profile-likelihood-based fesl slafislic [52]. 
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Figure 1: The missing transverse momentum (E™®®) distributions of the 2-jet events in the signal region for the 
(a) Q-b-tag, (b) \-b-X.ag and (c) 2-Z7-tag categories. The data are compared with the background model after the 
likelihood fit. The bottom plots show the ratio of the data to the total background. The signal expectation for 
mn = 125 GeV and BR(H inv.) = 100% is shown on top of the background and additionally as an overlay line, 
scaled by the factor indicated in the legend. The total background before the fit is shown as a dashed line. The 
hatched bands represent the total uncertainty on the background. 
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Figure 2: The missing transverse momentum (E™®®) distributions of the 3-jet events in the signal region for the 
(a) Q-b-tag, (b) \-b-lag and (c) 2-b-iag categories. The data are compared with the background model after the 
likelihood fit. The bottom plots show the ratio of the data to the total background. The signal expectation for 
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Figure 4: The dijet invariant mass (nihj) distributions in the signal region for the l-Zi-tag category, for events with 
in the range (a) [120-160 GeV], (b) [160-200 GeV], (c) [200-300 GeV] and (d) [> 300 GeV], The data are 
compared with the background model after the likelihood fit. The bottom plots show the ratio of the data to the total 
background. The signal expectation for mn - 125 GeV is shown on top of the background and additionally as an 
overlay line, scaled by the factor indicated in the legend. The total background before the fit is shown as a dashed 
line. The hatched bands represent the total uncertainty on the background. 
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Figure 5: The dijet invariant mass (mth) distributions in the signal region for the 2-fo-tag category, for events with 
in the range (a) [120-160 GeV], (b) [160-200 GeV], (c) [200-300 GeV] and (d) [> 300 GeV], The data are 
compared with the background model after the likelihood fit. The bottom plots show the ratio of the data to the total 
background. The signal expectation for mn - 125 GeV is shown on top of the background and additionally as an 
overlay line, scaled by the factor indicated in the legend. The total background before the fit is shown as a dashed 
line. The hatched bands represent the total uncertainty on the background. 
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Figure 6: Upper limits on (Tvh x BR(// —»inv.) at 95% CL for a Higgs boson with 115 < mn < 300 GeV. The full 
and dashed lines show the observed and expected limits, respectively. 

At niH = 125 GeV, for VH production, a limit of 1.1 pb is observed compared with 1.1 pb expected. 
These combined results for VH production assume the SM proportions of the WH and ZH contributions. 
Observed (expected) limits are also derived for the two contributions separately, 1.2 (1.3) pb for WH and 
0.72 (0.59) pb for ZH. As shown in Table 4, the 2-tag categories are almost only sensitive to ZH, the 
1-tag categories are equally sensitive to WH and ZH, and the 0-tag categories are more sensitive to WH 
production. The two processes contribute approximately equally to the sensitivity. 

For the discovered Higgs boson at tue - 125 GeV, an observed (expected) upper limit of 78% (86%) 
at 95% CL on the branching ratio of the Higgs boson to invisible particles is set. These limits are de¬ 
rived assuming SM production and combining contributions from VH and gluon-fusion processes. The 
gluon-fusion production process contributes about 39% (29%) to the observed (expected) combined sens¬ 
itivity. 


8 Summary 


In summary, Higgs boson decays to particles that are invisible to the ATLAS detector are searched for 
in the final states of two or three jets and large missing transverse momentum m & pp collision dataset 
corresponding to an integrated luminosity of 20.3 fb“^ at a centre-of-mass energy of 8 TeV No excess 
of events over the expected backgrounds is observed. The results are used to constrain the cross section 
for VH production followed by the decay H inv. for 115 < mn < 300 GeV. The observed 95% CL 
upper limit on <tvh x BR(// ^ inv.) varies from 1.6 pb at 115 GeV to 0.13 pb at 300 GeV. Assuming SM 
production and including the gg ^ H contribution, an observed (expected) upper limit of 78% (86%) on 
BR(// ^ inv.) is derived for the discovered Higgs boson with tue = 125 GeV. This independent result is 
comparable to that of the ATLAS ZH search with Z ^ (i and H inv. [19]. 
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